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ABSTRACT 

We present a general procedure for obtaining the present density fluctuation 
probability distribution given the statistics of the initial conditions. The main 
difficulties faced with regard to this problem are those related to the non-linear 
evolution of the density fluctuations and those posed by the fact that the fields 
we are interested in are the result of filtering an underlying field with structure 
down to scales much smaller than that of filtering. The solution to the latter 
problem is discussed here in detail and the solution to the former is taken from 
a previous work. 

We have checked the procedure for values of the rms density fluctuation 
as large as 3/2 and several power spectra and found that it leads to results 
in excellent agreement with those obtained in numerical simulations. We also 
recover all available exact results from perturbation theory. 



Subject headings: 
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1. Introduction 

The current model for the formation of large-scale structure is based on the gravitational 
instability of small initial fluctuations in the energy density. These are assumed to have 
originated from quantum fluctuations during the inflationary epoch. The fluctuations form 
a random field, which is fully specified by the n-point joint probability density functions 
(hereafter PDFs). The one-point PDF is a useful tool for testing the statistical character of 
the initial conditions (i.e. whether they are Gaussian or others), since its evolution does not 
depend on the nature of the dark matter (Trimble 1987). So it would be rather convenient 
to develop an efficient and accurate procedure for obtaining the present PDFs given that 
for the initial conditions. Many attempts have been made at solving this problem. Hoffman 
(1987) used the Zel'dovich (1970) approximation for computing the variance of the density 
fluctuations in the mildly non- linear regime (a 2 < 1). However, he did not take into account 
the fact that evolution and filtering are non-commuting operations. To obtain the variance 
on comoving scale r, one should evolve the whole field (with all its structure) and then filter 
it on scale r. But in the analytical procedure followed by Hoffman the field is first filtered 
and then evolved. In the linear regime (a C 1) both procedures lead to the same result. 
However, for a ~ 1 the results obtained by these procedures are rather different. 

A procedure for obtaining the present PDF given the statistics of the initial conditions, 
that accounts for non-commutability of evolution and filtering was first described by one of 
us (Betancort-Rijo 1991). This procedure was the result of taking into account the simple 
fact that overdense regions with present comoving size r correspond to larger comoving 
regions in the initial field. However, in developing this idea a couple of subtle points were 
involved that were not properly dealt within that paper. 

During more than a decade that has passed since these first attempts at dealing with 
the problem under consideration (the article by Betancort-Rijo was in fact written in 1988) 
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much attention have been given to it. New methods have been developed for approximating 
the non-linear evolution (Buchert 1989; Bouchet et al. 1995) and for the computation 
of statistics of the evolved field. Bernardeau (1994a) developed a procedure for taking 
into account the effect of filtering on the moments of the PDF in the small variance limit 
and used it to compute the first few hierarchical amplitudes. Then, assuming that these 
quantities do not depend on a, he obtained (Bernardeau & Kofman 1995) the present PDF 
for any value of a. However, this assumption, although it holds for some initial conditions, 
is not universal (as we shall see) and, even when it holds, it is not easy to justify a priori. 

Despite all the attention that has been devoted to the present problem, the idea 
described by Betancort-Rijo (1991), which is at the root of the possibility of developing a 
simple and direct (it is not necessary to obtain the vertex generation function) procedure 
for computing the present PDF in the non-linear regime analytically, has not been followed. 
The exception is an article by Padmanabhan & Subramanian (1993), but there the problem 
is treated in the same way as was by Betancort-Rijo (1991), thereby having the same 
drawbacks. 

In this paper we shall show how to correct the errors included in previous work and 
give the correct procedure for obtaining the present PDF. To check the correctness of the 
procedure we compute the values of S3 and £4 in the low-a limit (that have been computed 
analytically by Bernardeau 1994a) for several power spectra. Our procedure should render 
the exact values of all hierarchical amplitudes in this limit, and this is what is found for 
these two amplitudes. We also study the dependence on a of these amplitudes and of the 
variance in the low limit and for finite values of a. We then compute the PDF explicitly for 
several cases for which numerical simulations are available and find them to agree with the 
numerical results. 

In §3 we briefly review the basic idea (fully described in Betancort-Rijo 1991) involved 
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in the present approximation and show how to use it to derive correctly the present PDF 
given the initial one. In the second section we shall discuss the technique that we use for 
the treatment of the non-linear evolution. In the last section we present and discuss the 
results of several explicit calculations. 

2. Non-linear evolution 

The complete Zel'dovich approximation (CZA; Betancort-Rijo & Lopez-Corredoira 
2000) is an approximation to the non-linear evolution of density fluctuations depending 
only on three quantities, Aj, defined so that that 1 — Aj are the eigenvalues of the tensor ||, 
where u stands for the peculiar velocity of matter as given by the linear theory in comoving 
units, i.e. the peculiar velocity v = R(t)u; R(t) being the scale factor of the Universe 
and q stands for the Lagrangian coordinates, which are equal to the initial comoving 
coordinates (Eulerian). This tensor we call the linear local deformation tensor. The full 
local deformation tensor is defined in the same manner but using the exact value of u. This 
approximation is exact to the second order. An approximation exact to the fourth order 
may be obtained out of it but this implies introducing additional variables. It takes the 
form (Betancort-Rijo & Lopez-Corredoira 2000) 

(1 + 5)- 1 = n(l - ri(X)Xi + fi(X)xi + f 2 (X) Ul ), (1) 

i 

where fi, f'2, are certain functions of the variables Aj. The variables x, u which are 
given at any point by certain integrals over the whole field (Betancort-Rijo 2001) but, 
for statistical purposes we only need their probability distribution for fixed values of the 
Aj. In fact, to obtain the moments exactly to the third order and almost exactly to the 
fourth we can use expression ([!]) with u equal to zero and x equal to its average value for 
fixed values of the Aj (< Xi > (A)). This is the approximation that in principle we use. It 
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has a dependence on the power spectra through the mean value of x (see Betancort-Rijo 
& Lopez-Corredoira 2000). However, the dependence is so small that may be neglected 
without practical lost of accuracy. So, what we do in practice is to use simply the CZA 
(expression (1) with x — u — ), which is strictly exact to the third order for a particular 
power spectra (such that < X{ > (A) = 0) and practically exact for all of them. The CZA 
gives the density fluctuation at a point in term of the values of the Aj at that point. When 
dealing with filtered fields, however, the relevant quantity is the mean density within certain 
region (of any prescribed shape). It may be shown that this mean density fluctuation is also 
given exactly at least to the third order by expression ([|) with the Aj representing now the 
linear eigenvalues of the deformation tensor corresponding to the filtered field; that is, the 
deformation tensor derived from the peculiar velocity field corresponding to the growing 
mode of the filtered field. The x variables are now defined in a different manner. For 
non-filtered fields the presence of this variables is due to the different rate of growth of the 
various contributions to the local tidal field. In the present case, however, these variables 
have an extra component due to the filtering process. But, despite of the different origin of 
the x variables, expression ([I]) and the probability distributions for the x retain their form 
at least up to third order. So, setting in this expression x equal to its mean value for fixed 
Aj leads also to an approximation exact to third order (for statistical purposes) for filtered 
fields. There is, however, an important difference between the two mentioned contributions 
to the x, because the part due to the tidal field remains small for any power spectra while 
that due to the averaging process grows indefinitely with increasing power on small scale. 
This implies, as we shall find, that for power spectra with too much power on small scale 
the procedure presented here fails even for values of a smaller than one. However, even 
in these cases the results are exact to third order, but many orders contribute for finite 
but small a . So, the steeper the power spectra the smaller the value of o at which the 
approximation start failing. This happens, however, only for very steep power spectra not 
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very relevant to cosmology. To deal with these cases we should use expression ([I]) with the 
mean values of x, not simply the CZA. 

3. PDF of the filtered field 

The effect of filtering on the field of density fluctuations have been treated exactly by 
Bernardeau (1994). In this treatment one may, in principle, obtaining (previously to caustic 
formation) all moments of the PDF of the field of density fluctuation filtered in some scale 
as a series of powers of a and sometimes (when the corresponding Hedgeworth expansion 
converges) to reconstruct the PDF. However, even when this works it is complicated. The 
procedure described here to deal with this problem ,which is basically that in Betancort-Rijo 
(1991) is more immediate, allowing the direct evaluation of the PDF even for values of a 
well above unity. Furthermore, it provides a clearer visualization of the effect of filtering. 
We now describe the basic idea of our procedure using the CZA, although the procedure in 
itself is independent of the approximation used for the evolution. 

Consider certain geometrical body, B\, centered at some point within the initial field. 
By assigning to this point the mean value of the density fluctuation within B\ and repeating 
the process for each point, we generate a new field which we call the field filtered within B\. 
We may now obtain the linear peculiar velocity field corresponding to the growing mode of 
this field. Let Xi(B\) be the proper values of the corresponding linear (linearly calculated) 
deformation tensor, and let B-i be the body which is the result of transforming B\ with this 
deformation tensor. The CZA, in the present context, must be written in the form: 

l + *(B 8 ) = (f[(l-r < (X(5 1 ))A i (S 1 )))- 1 . (2) 

i=l 

where 5 is the present value of the density fluctuations filtered with B 2 at the point under 



- 8- 



consideration. In our problem, the filtering body B 2 is a sphere of radius r and B\ is an 



It must be noted that in the actual evolution the present sphere does not transform 
(when evolved back) exactly into an ellipsoid. The real surface has the same quadrupolar 
component as the mentioned ellipsoid, but contains smaller wiggles. However, for the 
present problem those wiggles are completely irrelevant. 



approximation. By one-point expression we mean that the value of 5 at a given point (the 
mean value within a B 2 centered at that point) is expressed only in terms of the Aj at that 
point (the mean value within B\). To go beyond this approximation, we need to obtain the 
values of the \ at different points and integrate over them, or, what is the same, to include 
the x variables described in the previous section. 

It should be noted that in equation (|2|) the Aj depend on Bi, which, in turn, depends 
on the Aj (for a given B 2 ). It is possible to show that, as long as caustics do not form, for 
every point there is a unique ellipsoid (with certain axis and orientation) such that the local 
deformation tensor of the field filtered with the ellipsoid transforms that same ellipsoid into 
a sphere of radius r (B 2 in general). This prescription uniquely determines the values of the 
Aj (hence, 5) at any given point. These unique Aj (for a given B 2 ) we call Aj. We must now 
obtain the PDF for the Aj at a randomly chosen point in the initial field. With this PDF 
and equation (Q), the derivation of the PDF for S is a standard problem. 

For a field filtered with a given body (the same at every point; usually a sphere of given 
radius), the PDF for the Aj is given by (Doroshkevich 1970): 



ellipsoid with semiaxis r(l — Tj(A)Aj) 



. Hereafter, we write simply Aj instead of Xi(Bi). 




expression for 5 as given by the CZA 



P(X 1 ,X 2 ,X 3 ,a) 



5 5 / 2 27 



(A 1 -A 2 )(A 1 -A 3 )(A 2 -A 3 ) 



87TO- 1 



•6 




) 
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o 2 = (5 2 L ) = ((Ai + A 2 + A 3 ) 2 ); A 3 < A 2 < A x . (3) 

a 2 is the variance of the initial field of density fluctuations, Sl, filtered with the 
given body and extrapolated linearly to the present time. The trace of the deformation 
tensor extrapolated linearly to the present time (is given the tensor used in the Zel'dovich 
approximation) by 5l- Betancort-Rijo (1991) and Padmanabhan & Subramanian (1993) 
took the PDF for the \[ to be simply proportional to the PDF for the Aj (with Aj = AQ for 
the field filtered with the ellipsoid that will contract to a sphere of radius r; that is, with 
an ellipsoid with semiaxis r(l — rj(A')A^) _1 . Note that now the filtering body changes from 
one point to another, so a is no longer a constant but a function of the A^. 

Ax.r) = (^55/ W{FfP(K)d 3 K- 
= 3[sinF — FcosF] 

pa 

with F = ^/J2i K 2 E 2 r- E { = (1 - ^(A')^) -1 , K = \K\; P(K) = power spectra. 

In practice, in the papers mentioned and, to some extent, here, instead of the full 
equation @ a more convenient approximation is used. Since a depends mainly on the 
volume of the filtering body (this is exactly so for white- noise power spectra), only the 
dependence on this quantity is considered. We then have: 



a 2 (X\r) = ^ 2 (0>(n(l-r 4 (A~OAr 1/3 )) = ° 2 (r(l[(l- n (A') A'J ~ 1/3 ) ) = a 2 (r(l + 5)^), (5) 
i=i i=i 

where a 2 (r) is the variance of the linear field filtered with a sphere of radius r. However, if 

we want to obtain results which are exact to the fourth order, as we shall do in some cases, 

the small shape dependence of <r(A,r) must be taken into account. In these cases, we use 

the following excellent approximation to the actual shape dependence of a: 
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cr 2 (A, r) = cr 2 (r(l + 5) 1 ^ 3 )Z 2 (X) 



Z(X) ~ exp 



(6) 



To compute the constant B in the exponential we may use the exact form of Z in the 
case when two of the A are equal. In this case and for power law power spectra (with index 
n) expression (R) reduces to: 



Z(P,n) 



P'^ Jo°°io°> 2 + P 2 s 2 )%(W(y/v? + s 2 )) 2 s ds du 



$™W 2 {u)v?+ n du 



P 



1- A 3 
1-A X 



Ai — A2 7^ A3 



(7) 



We find: B(n = -2)=0.0486; B(n = -1)=0.0499; B(n = 0)=0. 



The approximation used in the mentioned works (which use expression (5)) may be 
expressed in the form: 



P'(A>) = C P(A,a(r(l + 5(A)) 1 /3))| x=A ,, (g) 

where P' is the PDF of the A£; 5(A) is given by equation (fj) (with x=u=o); a stands for 
cr(r), and C is a normalization constant. This constant has to be introduced because, unlike 
P, the integral of P' over all A' space is not equal to unity, due to the dependence of a on A. 

We now show why equation is incorrect and give the correct one. To this end 
we consider an approximation in which 5 depends only on 5l(= X\ + A 2 + A3). We do 
this because the relevant fact at this point is basically the same with one and with three 
variables, but it is much easier to understand in the former case. 

For definiteness, let us assume that: 
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) 



-3 



1 + 5 



(9) 



an approximation that has been used in Betancort-Rijo (1991). 

For Gaussian initial conditions, Sl, which corresponds to a fixed filtering scale, follows 
a normal distribution. To obtain the PDF for 5' L (the value of 5l within the sphere of 



given point 5' L takes the value 5° L means that at this point the linear value of 5 filtered with 
a sphere of radius rj (1 — 5^/3) is just 5° L . The accumulated probability, P\5' L > 51), that 
at a point chosen at random in the initial field 5' L is larger than or equal to some given value 
may be obtained immediately. To this end we only need noting that if a sphere of radius 
r/(l — 5' L /3) (with 5' L > 51) have to contract to a sphere of radius r and shell crossing does 
not happen, the sphere of radius r/(l — 51/3) will have to contract to a radius smaller than 
r. So, the value of 5l within the sphere with radius R/{1 — 51/ 3) must be larger than 51, 
since if the value of 5l within this sphere were just 5° L it would contract to a sphere of radius 
r. Thus, the accumulated probability under consideration must be equal to the probability 
that the value of 5l within a sphere of radius r/(l — 5° L /3) be larger than or equal to 51- 
So, we have 



where erfc stands for the error complementary function. To obtain the probability 
distribution for the values of 5' L (the value of 5l within the sphere that will contract to 
a sphere of radius r) at a randomly chosen point we only need deriving the cumulative 
probability with respect to 5° L . 



radius r/(l — %) that have contracted at present to the sphere of radius r) that, by means 
of equation (||), allows us to obtain the PDF for 5, we proceed as follows. To say that at a 





(10) 
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P'(S'l 



d 



d5° r 



P'(h>5° L )\ 



P'{5' L )d5' L 



1 „-62/(2* 2 (6' L )) d ( 6 L 



(11) 



Note that since a depends on 5' L we cannot take it out of the differential. So we cannot 
write: 



P'(S' L ) = CP(6 L ,a(6 L ))\ SL=s , L 

contrary to the assumption in equation (|§p. With these considerations in mind, it is easy to 
obtain the correct distribution for the A-. We need only to generalize the above argument 
to three dimensions. The key argument is that the variables A^/cr(A') must follow the 
distribution followed by the variable Aj/cr for fixed filtering (when a is a constant). So we 
have: 



P'(A>)= P(A>(A)V 3 (A) 



a(X)=a(X,r). 



U=A" 



(12) 



The last factor in this expression is the Jacobian of the transformation from the 
variables Aj/<r(Aj) to the variables A«. Note that the P' obtained by means of this expression 
is automatically normalized. 



Equation (12) is exact for any initial conditions (provided that the A, determine the 
evolution), although we shall use it only for Gaussian initial conditions; that is, with P 
given by equation (H). 
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P' is the PDF of the A^ (the Aj associated with the unique ellipsoid that will contract to 
the sphere of radius r centered at the given point) for a point chosen at random in the initial 
field; that is, in Lagrangian coordinates. However, we want to obtain the PDF for 5 at a 
point chosen at random in the present "physical" coordinates (i.e. Eulerian coordinates). 
So, we must first obtain the PDF for the A£ in Eulerian coordinates. To this end, we only 
need to multiply P' by the Jacobian of the transformation from Lagrangian to Eulerian 
coordinates. For a non-filtered field this Jacobian is simply 



n(l-r 4 (X)Ai). (13) 

i=l 

However, for filtered fields this expression is not correct. This can be appreciated 
by noting that the resulting distribution for the A- (the result of multiplying P' by the 
Jacobian in question) is not normalized. For non-filtered fields, on the other hand, it is easy 
to check that the Jacobian given by equation fll3|) leads to a normalized PDF. The correct 
Jacobian must exhibit the following properties: it has to be zero at points where 5 goes to 



infinity; it must reduce to expression fll3|) in the limit of small power on scales below r; it 
must render a normalized PDF with vanishing first moment for any power spectra. The 
following expression meets these requirements and is exact at least to first order 



^(A) = n( 1 -( 1 -7h(A)A J ), (14) 

i 

where 7 is certain spectral constant that tends to zero in the limit of small power on 
scales below r. We cannot say a priori whether this expression is exact to higher orders. 
However,it seems to be exact to second order (within the precision of our computations) 
and it gives excellent results to all orders. It must be noted that expression (O) is not the 
value of the Jacobian at a point with the given Aj. This would be so if the evolution were 
determined solely by these quantities. However, we know that this it is not strictly true, so 
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expression (|14]) should represent the mean value of the Jacobian over all points with a given 
value of the Aj. To obtain 7 we have considered a fluctuation with spherical symmetry 
and the evolution of the region such that the mean density fluctuation within a sphere 
centered at points belonging to it is above some threshold. Using the spherical model with 
the mean density profile, we obtained the contraction of this region and by comparison 
with expression (|T4[) we find 7. However, this derivation is lengthy and although it might 
be interesting in itself, it is not essential to the present purposes, since we may obtain that 
constant using any relevant analytical result available (i.e., the leading order value of S3) 
and comparing it with the one obtained using expression ([[4]). In either case, we find: 



l d[ln(a\r))\ 

1 6 d[ln(r)] y ' 



4. Results and conclusions 

To check the accuracy of the approximation developed in this article, we shall present 
some explicit calculations and compare them with the results of the corresponding available 
numerical simulation and with exact results obtained with the perturbation theory. We 
shall consider Gaussian initial conditions with power spectra of the form P(k) = Ck n . 



In these cases, equation (EH takes the form: 



a(X,r)=a (l+5)' a / 2 ; 

(T = a(r);a = , (16) 



with 5 given by the CZA 



(1 + = na - 

i 



(17) 
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To derive exact results we use for rj(A) its exact expression (Betancort-Rijo & 
Lopez-Corredoira 2000). To derive result for finite values, however, we find it much more 
convenient to use the following compact approximation 



g(x,n) 



r i (A) = </(?(A J + A jr ),n) 



o / n2/63 \ 0.5533 

! (i - *?sr) V < 1-584 



1.55 



■x > 1.584 



(18) 



Using this in expression (12), with P(X, a) given by (3) and multiplying it by the 
Eulerian-Lagrangian Jacobian, J(A, a), (expression (14)), we obtain the PDF for the A^ 
corresponding to a point chosen at random in Eulerian coordinates, which we represent by 
G(AVo): 



G(A> )= P(A,a(A)V 3 (A) 







)J 


d\ 3 





J(X',a). 



(19) 



The moments of the PDF of the present value of the density fluctuation top-hat filtered 
on scale r are then given by: 



d\[ dX'J dX' 3 G(X',a )5 n , (20) 

-oo <j — oo J— oo 

where 5 is given by equations (17) with A = A', provided that each of the factors in (17) is 
positive and 5 is smaller than 170. If this expression leads to a value of 5 larger than this, 
we use the following expression: 



5 = 170 



1.589 



This last prescription represents an attempt to deal with caustic formation in a simple 



-16- 



manner, so that equation (23) gives sensible results for any value of a. However, it is not 
relevant for the results given here. 

The PDF for S, P(S, a ), may be obtained in the very process of computing any of the 
moments. To this end, we add up all the values of G(\, a) (multiplied by the corresponding 
volume of the integration step) that are found during the integration to lie between 5 and 
S + AS (for a suitably chosen AS, depending on S) and divide the result by AS. 

The values of S given by (17) for very under dense regions is not correct this imply 
(S) ^ 0. However, this problem disappears if we set P(S) = for S < —0.82. 

The values of S3 and £4 given by 

a _ (^) . o = (fl 4 )-3(J 2 ) 2 

3 (<5 2 ) 2 ' 4 {S 2 } 3 

in the low-a limit, may be computed analytically. The exact value of 
S 3 (a = 0) = (34/7 — (n + 3)) is obtained, for the power spectra under consideration, when 
expression (17) for S is used; the value 4 — (n + 3) is found when the are set equal to 
one. For a non-filtered field we have computed S4 (a = 0) analytically both for the CZA 
(equation (17)) and for the Zel'dovich approximation. Note that the non-filtered case may 
be treated formally as a filtered field with negligible power below the filtering scale, so that 
we may use expression (30) with a independent of A and 7 = 0. For power law power 
spectra this corresponds to the limit n = —3, but it corresponds more generally to any 
power spectra having the mentioned properties. The correct values (60712/1323; 279/9; 
Bernardeau 1994a) are found. For the evolution of these amplitudes and the variance in the 
small a limit we find 



(S 2 ) = a 2 + 1.81a 4 + 3.8a 6 
S 3 = 4.86 + 9.815a 2 + 48a 4 
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S 4 = 45.89 + 288.8a 2 + 2072a 4 



(21) 



This results are exact to next leading order for S3 (in principle only for n = —2.7 
but the n dependence is negligible) and almost exact to the same order for (S 2 ), £4. As 
pointed out in section 2, we could use the statistical properties of the x, u variables to 
obtain these two quantities exactly to the stated order, but the second coefficient in the 
expansion for (5 2 ) will increase by less than a per cent which is about the numerical error 
for this coefficient. So, it would be an unnecessary complication to use here those additional 
variables. Furthermore, the purpose of the present procedure is to obtain directly and 
accurately the PDF for finite values of a. The question of the maximum order of exactness, 
although interesting, is only of peripheral relevance to our purposes. These results agree 
with those exact results given by the perturbation theory by Scoccimarro & Frieman (1995) 
and the remaining ones are close to those obtained by Fosalba and Gaztanaga (1998) using 
the spherical model. We see that this model predict well the hierarchical amplitudes but 
not so well the moments. The above results provide a check to the CZA but not properly to 
the procedure presented here, since for non-filtered fields expressions (12) and (14) reduce 
to the standard form. A check to these expressions is provided by the cases of power law 
power spectra with n > —3. For n = — 2 we find: 



(S 2 ) =a 2 + 0.88a 4 + 0.54a 6 



S 3 = 



3.86 + 3.12a 2 



S 4 = 27.56 + 58a 2 



(22) 



Here, to obtain the second coefficient in (5 2 ) exactly we have considered the small 
(=0.015) contribution due to the dispersion of the x variables. This is not necessary for 
the amplitudes. It is also important to consider the shape dependence of a (see expression 
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(6)), otherwise we would obtain a coefficient 0.1 too large. For n = —1, n = 1 we find 
respectively: 

(S 2 ) = a 2 + 0.223a 4 

5 3 = 2.86 + 0.267a 2 

5 4 = 13.89 - 2.583a 2 

(S 2 ) = a 2 -0.173a 4 

S 3 = 0.86 + 0.8a 2 
S 4 = 0.556 - 17.8a 2 (23) 

For n > —1.5, S 3 and S4 are practically independent of a. For example, for n = —1.3 
the change of these quantities up to a = 1.52 is less than a 4%. For smaller values of n, 
however, the evolution of this quantities is more apparent. For n = —2 we find (table 1) that 
S3 and £4 increase considerably from a = to a pa 0.5 and more steeply for higher values 
of a. In figure 1 we present the PDF for 8 corresponding to P(k) oc k" 1 and a = 0.3; 0.5. 
Comparing this with figure 7 of Bernardeau & Koflman (1995), we find them to be almost 
identical. The PDFs of Bernadeau & Koffman were obtained under the assumption that 
the S p parameters (the p-th cumulants divided by a 2 ( p ~ 1 ' ) ) are independent of a. We have 
seen that this is the case for n = —1, at least for S 3 and S4. So, the agreement between 
both pairs of PDFs is to be expected. The evolution of {S 2 ) agrees very well with the result 
by Lokas et al. (1996). 

In figure 2 the PDFs corresponding to a = 0.92; 1.52, n = —1.3 are represented. These 
should be compared with figure 10 of Bernardeau & Kofman (1995). The agreement is 
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very good, our results being always within the numerical errors of the results given in that 
article. 

It should be noted that the mentioned numerical results assume a cold dark matter 
(CDM) while we have used a power-law power spectrum with the same local logarithmic 
derivative. Using the whole CDM power spectrum in our computations would represent 
only a small complication. If we have not done so, it is because this implies only a slight 
improvement in the results that is completely unobservable given the precision of the 
simulations. 

The procedure described here for obtaining the PDF for 5 is in principle valid only 
for top-hat filtered fields. However, the basic lines of reasoning followed here are entirely 
translatable to another smoothing procedures. 

For a given smoothing (in the present field), we must determine the smoothing in 
the initial field rendering the minimum dispersion for the smoothed present field. That 
is, we must find a smoothing procedure that, when applied to the initial field at a given 
point, renders a set of A- (extrapolated linearly to the present time) that predict the mean 
value of 5 at the transformed point, given by an expression similar to (17), and so that 
the dispersion in 5 for given A-s is a minimum. This is what we have implicitly done here 
for the spherical top-hat smoothing, and we have found that it corresponds to an initial 
ellipsoidal top-hat smoothing. 

The solution to this problem for Gaussian smoothing is not simple, although one could 
guess that a triaxial Gaussian, with the axes determined as in the top-hat case, should not 
be a bad approximation. Within this simple approximation and for a power-law spectrum 
the PDF for 5 must be the same as for the top-hat case with the same a. For more general 
power spectra, we cannot decide a priori whether the value of the logarithmic derivative 
must be taken on the scale where the top-hat has the same volume as the Gaussian or on 
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the scale where it has the same a. 

As an example we discuss one of the cases represented in figure 4 of Kofman et al. 
(1994), where the Gaussian smoothing radius is R s = 6/i~ 1 Mpc; a = 0.55. If we matched 
the volumes, the corresponding top-hat radius is 9.38 /i _1 Mpc, so that for a CDM power 
spectrum n pd —1. The PDF corresponding to a top-hat filtering with a = 0.55; n = —1 is 
represented in figure 3 and shows very good agreement with the numerical results in the 
last reference. If we matched the cr's, the corresponding top-hat smoothing scale imply 
n = —1.1. The corresponding PDF is almost indistinguishable from the former. 

We have seen that to obtain the PDF for 8 we must integrate over all values of A 
subject to the condition 8(X) = 8. If we considered, however, an approximation where 8 
depended on the Aj only through 8l = X\ + A2 + A3, the computation of P(8) would be 
greatly simplified. We have already seen that such an approximation cannot be entirely 
correct, but leads to a simple analytical expression that displays, in a clear manner, the 
effect of smoothing. Taking for 5(Sl) that corresponding to spherical collapse and using the 
fact that the distribution for 8l is simply a Gaussian, equation ([11]) and (14) immediately 
gives for P(S): 



P(S) 



((l + 5)- a602 + f^(5)) 




a = 



n + 3 
3 



6 L (6) = 1.66(1- (1 + 8) 



-0.6024 



(24) 



Where the last equation represents a very accurate approximation to the spherical 
model that we propose. However, expression (|24|) represents only a rough approximation to 
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P{5). For a = 1.52 the value of P(S) given by (|24"D may differ from the real value by a factor 
greater than three. Furthermore, the value of S 3 (a = 0) given by ([24]) is 37/7 — (n + 3), 
which is 3/7 too large. 

We may conclude that for 5 G (—0.9, 10) and a as large as 1.52, the P(5) computed 
as indicated here is, so far as we have been able to check, in excellent agreement with 
the results of numerical simulations. This agreement is particularly meaningful since our 
derivation of P(S) involves neither free parameters nor modeling freedom (except for 
some minute refinements non-essential for the results given here). Equation ([12]) is exact 
and is valid for any field. Equation (14) for the Jacobian of the Eulerian-Lagrangian 
transformation in a filtered field is also exact (at least to first order in Aj) and, although 
derived for a Gaussian field, it seems to be valid for any field. 



-22 - 
REFERENCES 

Bernardeau, F. 1994a, ApJ, 433, 1 

Bernardeau, F. 1994b, ApJ, 427, 51 

Bernardeau, F., & Kofman L. 1995, ApJ, 443, 479 

Betancort-Rijo, J. 1991, MNRAS, 251, 399 

Betancort-Rijo J., Lopez-Corredoira M., 2000, ApJ 534, L117 

Betancort-Rijo J., 2001, in preparation 

Bouchet, F. R., Colombi, S., Hivon, E., & Juskiewicz. R. 1995, A&A, 296, 575 

Buchert, T. 1989, A&A, 223, 9 

Doroshkevich, A. G. 1970, Astrofizica (VI) 6, 4, 581 

Fosalba P, Gaztahaga E., 1998, MNRAS 301, 503 

Hoffman, V. 1987, ApJ, 318, L7 

Kofman, L., Bertschinger, E., Gelb, J., Nusser, A., & Dekel, A. 1994, ApJ, 420, 44 
Padmanabhan, T., & Subramanian K. 1993, ApJ, 410, 482 
Scoccimarro R., Frieman J., 1995, ApJ 473, 620 

Tokas E. L., Juszkiewicz R., Bouchet F. R., Hivon E., 1996, ApJ 467, 1 
Trimble, V. 1987, ARA&A, 25, 425 
Zel'dovich, Y. B. 1970, A&A, 5, 20 

This manuscript was prepared with the AAS L^T^X macros v4.0. 



-23- 

FIGURE CAPTIONS 

Fig. 1: Probability distribution function (PDF) for o 
Fig. 2: Probability distribution function (PDF) for o 
Fig. 3: Probability distribution function (PDF) for a 
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